Dynamical correlations near dislocation jamming 
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Dislocation assemblies exhibit a jamming or yielding transition at a critical external shear stress 
value a — a c . Nevertheless the nature of this transition has not been ascertained. Here we study the 
heterogeneous and collective nature of dislocation dynamics within a crystal plasticity model close 
to cr c , by considering the first-passage properties of the dislocation dynamics. As the transition is 
approached in the moving phase, the first passage time distribution exhibits scaling, and a related 
peak dynamical susceptibility xl diverges as xl ~ (°~ ~ °~c)~ a , with a « 1.1. We relate this scaling 
to an avalanche description of the dynamics. While the static structural correlations are found to 
be independent of the external stress, we identify a diverging dynamical correlation length t; y in the 
direction perpendicular to the dislocation glide motion. 
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PACS numbers: 61.72.Lk, 68.35.Rh, 81.40.Lm 

The mechanical behavior of crystalline solids subject 
to an external shear stress a is controlled by the existence 
of a finite yield stress a c : for stresses below a c only elas- 
tic or reversible deformation of the material takes place, 
while sustained plastic or irreversible deformation medi- 
ated by dislocation motion is observed for a > a c . The 
basic phenomenology of two states with different rhe- 
ology is the same also for many amorphous materials 
ranging from foams to amorphous and granular media, 
even though the microscopic mechanisms responsible for 
a transition to "flow" are not yet well understood. The 
concept of jamming has been introduced to comprehend 
the observed phenomenology: In a finite region of the 
parameter space (spanned by control variables as tem- 
perature, density and an external force), due to a "cage 
effect", or self-induced constraints on the motion of the 
system constituent particles, the phase space dynamics 
gets restricted, i.e. jammed [l|-[3- I n particular for two- 
dimensional ensembles of crystal dislocations a yielding 
transition at a finite a c has been established, even in the 
absence of any external sources of disorder [5j. The cru- 
cial mechanism here is the generation of effective disor- 
der in the local stress field, due to constrained dislocation 
motion by the combined effect of the slip geometry and 
the long range anisotropic dislocation-dislocation inter- 
actions, which are also manifest in the formation of typi- 
cal, metastable elementary structures such as dislocation 
dipoles and walls. 

On approaching the jamming transition from the mov- 
ing phase, the dynamics generally becomes increas- 
ingly heterogeneous, and a growing dynamical correla- 
tion length characterizing correlation across intervals of 
both space and time has been observed for granular sys- 
tems 0, 0], as well as on approaching the glass transi- 
tion of molecular liquids and colloidal suspensions [8| . In 
this work, we consider dislocation jamming in a similar 
vein. Known features of the 2d discrete dislocation dy- 



namics (DDD) model studied here are a relaxation of 
the strain rate or an Andrade law close a = <r c , fol- 
lowing a power law e ~ t -6> , « 2/3 [5]. The steady 
state strain rate displays a non-linear rheological behav- 
ior, e ~ (<j — <j c )^, with f3 « 1.8 [5|. The model also re- 
produces the experimentally observed scale free size dis- 
tributions of avalanches of plastic deformation It is 
also similar to the DDD models used in the context of 
many engineering applications [lOj. The question is now, 
whether one can establish the jamming/yielding transi- 
tion as a second order phase transition. This would imply 
the existence of a length scale, the correlation length, 
diverging at the transition point. The dislocation sys- 
tem is particular, since here the transition is of a non- 
equilibrium kind: the temperature is included only indi- 
rectly via the dislocation mobility and is thus irrelevant. 

The avalanches exhibited by dislocation assemblies [9] 
are suitable candidates for fundamental "dynamical het- 
erogeneities" or localized events. We define appropri- 
ate statistical quantities to characterize the dynamics via 
time-dependent first passage probabilities which measure 
the likelihood that a dislocation becomes liberated from 
the confining stress field over an observation scale. The 
first passage time distributions are found to display scal- 
ing close to the transition, and analogously to other sys- 
tems exhibiting jamming, the peak value of a dynamical 
susceptibility and a dynamical correlation length grow on 
approach to jamming and diverge at a = a c . At the same 
time static correlations describing the dislocation struc- 
tures are found to be virtually independent of the exter- 
nal stress. We also connect the results to an avalanche 
description of the dynamics. 

The DDD model can be considered to represent a 2d 
cross section (xy plane) of a 3d single crystal with a single 
slip geometry [HI, 0, [ll[ . The dislocations are taken to be 
straight parallel edge dislocations, with the dislocation 
lines along the z axis. The N point-like dislocations glide 
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FIG. 1: (color online) Scaled first passage time distributions 
according to Eq. Q, for various values of a and /. The value 
of the exponent f3 is f3 = 1.8. 



along directions parallel to their Burgers vectors 6 n = 
s n bu x , with s n = ±1 and u x the unit vector in the x 
direction. Equal numbers of dislocations with positive 
and negative Burgers vectors are assumed. Dislocation 
climb is not considered for simplicity. The dislocations 
interact with each other through their anisotropic long- 
range stress fields, 
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where D = /x/27r(l — v\ with ji the shear modulus and v 
the Poisson ratio of the material. The dynamics is taken 
to be overdamped with a linear force-velocity relation, 
resulting in the equations of motion 
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where v n is the velocity of the nth dislocation, Xd is the 
dislocation mobility, s n is the sign of the Burgers vec- 
tor of the nth dislocation, and a is the applied external 
shear stress. The long range interaction forces are com- 
puted by imposing periodic boundary conditions in both x 
and y directions. The equations of motion are then inte- 
grated numerically with an adaptive step size fifth order 
Runge-Kutta algorithm. This is done in dimensionless 
units by measuring lengths in units of 6, times in units of 
l/(XdDb), and stresses in units of D. If the distance be- 
tween two dislocations of opposite sign gets smaller than 
26, they are removed from the system, a process reminis- 
cent of the annihilation reaction occurring in real plas- 
tically deforming crystals. The simulations are started 
from a random initial configuration of No dislocations 
within a square cell of linear size L. The system is first 
let to relax with a = 0, in the absence of external stresses, 
until it reaches a met ast able dislocation arrangement. 
During this initial relaxation, a significant fraction of the 
dislocations get annihilated. We consider here systems 
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FIG. 2: (color online) Top: The average self-overlap order 
parameter Q(l,r) for various values of a as a function r, for 
1 = 1. Bottom: The corresponding four-point dynamic sus- 
ceptibilities X4(Z,t). The inset shows evidence of a data col- 
lapse of X4(t), obtained by rescaling the data with the scaling 
forms of (tf p ) and X4*, with a = 1.1 and f3 = 1.8. 



of linear size L = 2006 and L = 3006, with 7V = 1600 
and 7V = 3600, respectively. After the relaxation, the 
external stress is switched on, and the evolution of the 
system is monitored. For a > a c ~ 0.025 the sys- 
tem eventually reaches a steady state, with the strain rate 
e ~ b n v n fluctuating around a constant mean value. 
This steady state consists of TV = 350 — 400 dislocations 
for L = 2006, and TV = 800 — 900 for L = 3006. 

To characterize the heterogeneous dynamics we first 
consider the first passage time distribution of dislocations 
in the moving steady state. The first passage time tf p 
is defined as the time at which a dislocation first moves 
across a distance I from some initial position, and thus 
becomes liberated from a confining stress field at this 
scale. Given the observed non-linear stress dependence 
of the steady state strain rate [5[, the mean first passage 
time is expected to behave like (tf p ) ~ ^ l(o~ — 

a c )~P . Assuming Stf p ~ (tfp), one can write 



P(t fp ]l,a) 



T 



tf p — I (a — a c ) 



. (3) 



Fig. [T] shows a data collapse according to Eq. (|3]) with 
f3 = 1.8, in good agreement with the scaling of strain 
rate with the applied stress [5| . The scaling function J-(x) 
exhibits a maximum and a tail which does not contribute 
much to the typical escape characteristics. 

A cumulative version of P(tf p ) has been called the in- 
stantaneous self-overlap order parameter [7| , 
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where N is the number of dislocations and w n equals 1 
if the displacement of dislocation n remains less than I 
across the time interval t —> t + r and equals otherwise. 
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The first and second moments 

Q(1,t) = (Q t (l,r)) 
Xi {l,t) = N[{Q t (l,r) 2 ) 
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of QtiJ-, t) (calculated from sample to sample fluctuations 
of Qt{l,r)) are then used to characterize the dynamics. 
As is customary for similar problems, we consider here 
only samples that are still active at the observation time 
[13|- By construction, Q(7, r) decays from one to zero as a 
function of r, while the four-point dynamic susceptibility 
X<±(h T ) exhibits a peak at an intermediate r = r|, and 
vanishes for both early and late times, see Fig. [2j Q* = 
Q(Z,r|) « 0.5 is nearly a constant. The peak value xt of 
X4(Z, t) depends on the proximity of the critical stress cr c , 
and has been argued to be related to the typical number 
of particles (here dislocations) in a correlated fast-moving 
domain or heterogeneity [7]. 

In the simulations, we observe a growing xt as G c is 
approached from above, with an apparent divergence at 
a = cr c ~ 0.025, see Fig. [3] For external stresses a < 
0.07, the stress dependence of xt can De characterized 
by a power law xt ~ ( a — °~c)~ a , with a ~ 1.1. This 
result does not exhibit any significant /-dependence. The 
scaling of xt can De connected to the avalanche dynamics 
close to the transition by directly computing the variance 
of Q(Z,r) at r|. Consider the distribution of avalanche 
sizes s with a <T-dependent cut-off, P(s) = s~ T f[s/(cr — 
<j c )~ 7 ]. By assuming that the number of dislocations n 
swept to a distance I during an avalanche obeys (n) ~ 
(s), and using the condition M(n) = N/2 (arising from 
the observation that Q* ~ 0.5) to estimate the number 
M of independent avalanches occurring up to the delay 



time r 4 *, we obtain xt = N(5Q*) 2 ~ N(^M{n)/N) 2 ~ 
(n)/2 ~ (cr — <r c ) -7 ( 2-r ). As (n) grows on approaching 
cr c , the number M of heterogeneities contributing to xt 
decreases: an estimate from the numerics for L = 3006 
leads to M « 19 for a 0.06 and M « 6 for a = 0.037. 
Given that the avalanche size exponent is known to have 
a value close to the mean field r~ 1.5 [11], a ~ 1 would 
indicate 7 ~ 2, again consistent with the mean field -like 
scenario 
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To further explore the possibility of identifying a di- 
vergent non-equilibrium correlation length on approach 
to dislocation jamming, we consider the conditional prob- 
ability gd(r,l) = P[w(l,T%,r) = 0\w(1,t£,0) = 0] that a 
dislocation at a position r has moved at least a distance 
I up to a delay time r|, given that a dislocation at the 
origin has done so. gd(r, I) provides a measure of the spa- 
tial structure of the typical heterogeneities contributing 
to xt- Figs. |4]a) and b) show examples of this dynam- 
ical correlation function for stress values a = 0.06 and 
a = 0.04, respectively. We focus here for simplicity on 
the case 1 = 1. Notice that closer to jamming, corre- 
lations clearly extend further in the y direction, while 
we could not identify a clear trend in the x direction. 
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FIG. 3: (color online) The peak value xt °f the dynamic 
susceptibility as a function of the distance from the critical 
stress for different system sizes L and length scale parameters 
/. The solid line corresponds to a power law xt 00 (cr — cr c )~ a , 
with a = 1.14. The inset shows the same data as a function 
of a with linear axis scale. 
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FIG. 4: (color online) The dynamic and static correlation 
functions from a system of size L = 3006 and 1 = 1, with the 
origin in the middle of the figures. In a) and b) we show the 
dynamic correlation functions gd(r,l = 1) for a = 0.06 and 
a = 0.04, respectively. In c) and d), the static correlation 
function g(r) for the same stress values (a = 0.06 in c) and 
a = 0.04 in d)) are displayed. 



At the same time, the static density-density correlation 
function g(r), computed from snapshots of the disloca- 
tion configurations in the steady state, does not show 
any significant dependence on a HE To estimate 

the dynamical correlation length, we fitted an exponen- 
tial function exp(—y/^ y ) to the ^(r)-data averaged over 
a narrow strip along the ^-direction. This leads to a 
divergent stress dependence of £ y ~ (a — a c )~ Uy , with 
v y « 1, see Fig. [5] A similar analysis in the x-direction 
is less conclusive: due to the tendency of dislocations 
to form wall-like structures, the statistics is significantly 
worse in the x direction, and no clear trend is visible 
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FIG. 5: (color online) Main figure: The dynamic correlation 
length ^ y in the ^/-direction, exhibiting a power law divergence 
close to cr c . Inset: The dynamic correlation length £ x in the 
x direction. The length scale parameter / = 1 is used. 



for the evolution of £ x with a (inset of Fig. [5j): £ x ap- 
pears to have a relatively large stress-independent value 
around ^ x « 356. Recall that in 3d DDD simulations [l7| . 
avalanches have been observed to be characterized by a 
lamellar shape, something that could be related to this 
large and stress-independent value of £ x . Notice also that 
the observed scaling of £ y is consistent with the mean field 
cut-off scaling of the avalanche size distribution, assum- 
ing that £ y measures the linear size of a typical avalanche 
in the y direction. 

In summary, we have shown that the dislocation jam- 
ming transition is accompanied by a divergent dynamic 
susceptibility and a divergent correlation length. This 
scaling, which we relate to the avalanche dynamics of 
the system, is purely dynamical in nature: the static 
correlation function characterizing the dislocation struc- 
tures does not exhibit any significant dependence on a. 
The observed exponent values a « v y « 1 indicate that 
the heterogeneities are compact and characterized by 
a single diverging correlation length in the y direction, 
thus clearly demonstrating that the dislocation jamming 
transition is indeed a true non-equilibrium second order 
phase transition. It would be interesting to test similar 
avalanche-based ideas also for plasticity of amorphous 
solids [181 . Notice that despite similarities with other 
relevant systems such as granular media, there are also 
significant differences: In the dislocation system the dy- 
namics stops due to the formation of various metastable 
dislocation structures (dipoles, walls, and more compli- 
cated multipolar structures) arising from the combina- 
tion of anisotropic long-range dislocation-dislocation in- 
teractions and restricted dislocation motion along a sin- 
gle slip direction. On the other hand, e.g. in granular 
systems jamming takes place when grains get stuck by 
short range contact forces within cages formed by their 
nearest neighbors, a process that is not accompanied by 
any large scale grain structures. The observation that in 
the present case the results seem to be independent of I 



close enough to the transition suggest that contrary to 
many other systems where I needs to be tuned to a spe- 
cific value related to an underlying "cage size", such a 
typical microscopic length does not exist here. This indi- 
cates that the same scale free property would apply also 
for the force landscape experienced by the dislocations, 
with the first passage properties related to the time it 
takes for the force landscape to deform on the observa- 
tion scale, due to the avalanche dynamics. 
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